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The subleading term of the heavy quark potential (the analogue of the Liischer term) is computed in a string 
model for the case of three quarks. It turns out to be positive in 2+1 dimensions, making the potential non- 
concave as a function of the scale for fixed geometry. The results are compared to numerical simulations of the 
lattice gauge theory. 



1. Motivation 

The potential of three heavy quarks has re- 
cently been the object of detailed numerical stud- 
ies |1I2) . which lend support to the so-called Y 
law, inspired by a string picture of heavy baryons. 
Here, we study the leading effect of string fluctua- 
tions on the potential, the analogue of the Liischer 
term |3I4I5| in the quark-antiquark potential. Un- 
derstanding such corrections to the Y law is par- 
ticularly important since the transition from A to 
Y law occurs at large quark separations £Q . 

2. Setup 

We study three static quarks in a D- 
dimensional Euclidean space with periodic Eu- 
clidean time extent T. The potential is obtained 
from the T — ► oo limit. 



^qqq = - 1™ 4 ln Z qqq 
1 — »oo 1 



(1) 



In the string picture, the classical ground state 
of the three-quark system is given by three strings 
meeting at a junction, whose position is deter- 
mined by the requirement of minimal total length 
of the strings. The balance of tensions implies an- 
gles of =§■ between the strings. The classical po- 
tential is proportional to the total length of the 
strings in this configuration, 

Ki = a(Li +L 2 + L 3 )+m = aL Y + m , (2) 

where a is the same string tension that appears 
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in the mesonic potential, m represents the self- 
energy of the junction. 

Following the mesonic case 0], we shall ex- 
pand the action to second order in the transverse 
fluctuations £ a (t, s) of the string world sheets 
a = 1,2,3 around the classical configuration. 
Transversality means = and e a ■ £ a = 0, 
where e a is a spatial unit vector in the direction 
of string a. With (j>(t) denoting the position of the 
fluctuating junction, the boundary conditions are 

L(t, 0) = x a , C(t, L a +e a -$(t)) = $ ±a {t) , (3) 

where (f>± a = 4> ~ e a (^a ■ </>), and periodic in t. 

Invariance under Euclidean transformations in- 
side the plane of the sheet and perpendicular to 
it fixes the leading term in a derivative expansion 
of the bulk string action to 



rfbulk 
J fluct 



E / d ^ ■ d ^ 

„_i JT a 



where T a is the domain of £ a implied by the 
boundary conditions ©. We also include a 
boundary term, 



Qbound 
D fluct 



dt\^ 



The change of area caused by fluctuations of the 
junction in the plane of a given sheet cancels in 
the sum over sheets. The three-quark potential 
including leading fluctuation effects can thus be 
extracted from the partition function 
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3. Calculation 

We shall compute Z qqq in two steps. For fixed 
</>, each of the string partition functions, 

Z a {$) = jB^ a e xp{-lJ\d£ a \ 2 } , 

splits into a minimal- area and a fluctuation part, 

Z a ($) = e-i ^ d ^^ 2 |det ra (-A)|-^ , 

where £ a ,min is harmonic and satisfies © and the 
determinant is computed with Dirichlet boundary 
conditions on the domain r a = {(f, s) | < s < 
L a +e a -<f(t)}. 

To leading order in 0, 



sinh(ajs) 
sinh(u;L a ) 



where (jj^ are the Fourier components of <p(t). 
This implies 

1 iCminl 2 = ]>>COth( W L a )|^| 2 +0(0 3 ) , 



which is the change in minimal area due to 4>j_ a - 
The determinant induces a renormalisation of 
both a and m. After applying a Pauli-Villars 
regularisation, the determinant, which can be ex- 
pressed in terms of the heat kernel of A, can 
be computed by mapping the domain T a confor- 
mally to a rectangle L' a x T. Since the conformal 
map cannot change the ratio L' a /T (the modu- 
lar parameter of the cylinder), one has to choose 
L' a = L a + f e a ■ 4>dt. To leading order in <p, 
the conformal map is 



f(z) = z + 



1 \ - e a ■ . 



VT ^ sinh(wL Q ) 



0(cj> 2 ) 



Following Ref. the determinant on T a can be 
related to that on the rectangle, and one finds 



lndetr (-A) = 2 In 77 
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where r\ is Dedekind's function. 

The Gaussian integral over <f> in Q can now be 
performed and V qqq be extracted from the limit 
T — > 00, cf. Q). One obtains: 



V — T/ rc 



r, 



l/L 



(D-$)V X ) L +0(L~ 2 ), 



V i/l ~ 24 
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Here, we have separated contributions from fluc- 
tuations in the plane of the three quarks, vj, L , 
and perpendicular to it, Vy L . Note that both 
expressions are homogeneous in L ai representing 
exactly the term. 

4. Check on the mesonic string 

As a check, we can (artificially) split the string 
connecting a quark and an antiquark into two sec- 
tions of length L\ and L2, pretending there is a 
junction in between. The mass of the junction 
should not affect the large-L behaviour, and we 
should recover the known result. In this case, 
there is no contribution Vy L and the integral in 
Vy L becomes 
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which just corrects the Liischer terms of the two 
string sections into that of the full string. 

5. Special cases 

The following special cases reveal interesting 
consequences of our result. In the equilateral 



3 



case, L a = L = it so happens that the fluc- 
tuations in the plane do not contribute. Those 
perpendicular to it yield 

Vi/l = V} /L + {p-Z)V$ L = -{D - 3) . 

This means that there is no l/L term in D = 3. 

Expanding about the equilateral case, L a = 
(1 + e a )L with J2 a £a = 0) one nn( ls 

a x a ' 

so in Z? — 3, the l/L term is positive. This makes 
Iqqq non-concave as a function of the scale L for 
fixed angles between the quarks. This is actually 
true for all geometries, not just almost equilat- 
eral ones. Fig. n shows a density plot of Vi/ L , as 
given by Eq. J5J, as a function of L\ and Li for 
fixed L3 = 1. In D > 4, the l/L term is always 
negative, so Vq qq is concave. 




Figure 1. ^V 1/L (L l ,L 2 ,l) in D = 3. 

6. Lattice gauge theory 

Numerical simulations of lattice gauge theory 
were performed in D = 3. Vq qq was extracted 
from Polyakov loop correlators using the method 
of Ref. [S] on a 48 2 x 32 lattice with lattice spacing 
a w 0.15 fm ((3 = 11), for one qqq geometry. The 
deviation from a pure Y law (with the measured 
mesonic string tension) is shown in Fig. [2] 

The asymptotic Y law is clearly approached 
from above as predicted, making the potential 
non-concave, but the magnitude of the l/L term 
appears too large: V 1 / L = 1.2(2) 24 ^ compared 
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Figure 2. Vq qq — a q qLy — c qq from lattice sim- 
ulations for rectangular, isosceles qqq geometries 
with either the large side or the small sides of the 
triangle aligned with lattice axes. The curves are 
fits c— b/Ly.The Y law is approached from above. 

to a coefficient of 0.39 as predicted by Eq. (JSJ. 
The discrepancy may be due to the small quark 
separations: even for the largest geometries, the 
shortest string in the classical configuration is 
only 0.45 fm long, so the string picture is hardly 
applicable. In addition, it should be noted that 
the numerical result is very sensitive to the es- 
timate of the mesonic string tension. Finally, a 
mandatory extrapolation to T — > 00 has not been 
performed yet, so that excited-state contributions 
may significantly affect the measured coefficient. 

A more thorough check of the string predictions 
should also include other geometries. Finally, it 
would be very interesting to extend the calcula- 
tion to larger SU(iV) gauge groups and compare 
with predictions from 1/N expansions. 
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